Raising Bi-O bands above the Fermi energy level of hole-doped Bi 2 Sr 2 CaCu 2 08+5 and 

other cuprate superconductors 
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The Fermi surface (FS) of Bi2Sr2CaCu20g+a (Bi2212) predicted by band theory displays Bi- 
related pockets around the (it, 0) point, which have never been observed experimentally. We show 
that when the effects of hole doping either by substituting Pb for Bi or by adding excess O in Bi2212 
are included, the Bi-0 bands are lifted above the Fermi energy (Ef) and the resulting first-principles 
FS is in remarkable accord with measurements. With decreasing hole-doping the Bi-O bands drop 
below Ef and the system self-dopes below a critical hole concentration. Computations on other Bi- 
as well as Tl- and Hg-based compounds indicate that lifting of the cation-derived band with hole 
doping is a general property of the electronic structures of the cuprates. 

PACS numbers: 74.72.Hs,74.25. Jb,71.18.+y,74.72.-h 



First-principles band theory computations on the 
cuprates have become a widely accepted tool for gain- 
ing insight into their electronic structures, spectral 
properties, Fermi surfaces (FS's), and as a starting 
point for constructing theoretical models for incorporat- 
ing strong correlation effects beyond the framework of 
the local-density approximation (LDA) underlying such 
calculationsi*2iii4. For example, in the double layer Bi- 
compound Bi 2 Sr 2 CaCu208 +( 5 (Bi2212) — perhaps the 
most widely investigated cuprate — the LDA generated 
band structure^ is commonly invoked to describe the 
doped metallic state of the system. Band theory how- 
ever clearly predicts the FS of Bi2212 to contain a FS 
pocket around the antinodal point M(tt,0) as a Bi-O 
band drops below the Fermi energy (Ef), but such FS 
pockets have never been observed experimentally 7 . This 
'Bi-O pocket problem' is quite pervasive and occurs in 
other Bi-compounds. 8 Similarly, Tl- and Hg-compounds 
display cation-derived FS pockets, presenting a funda- 
mental challenge for addressing on a first-principles basis 
issues related to the doping dependencies of the electronic 
structures of the cuprates. 

In this Letter, we show how the cation-derived band 
responsible for the aforemenentioned FS pockets is lifted 
above Ep when hole doping effects are properly included 
in the computations. Detailed results for the case of 
Bi2212 are presented, where hole doping is generated ei- 
ther by substituting Pb for Bi or by adding excess oxy- 
gen in the Bi-O planes. With 20% Pb doping in the 
orthorhombic crystal structure, the Bi-O band lies w 1 
eV above Ep and the remaining bonding and antibond- 
ing FS sheets are in remarkable accord with the angle- 
resolved photoemission (ARPES) measurements on an 
overdoped Bi2212 single crystal^. Below a critical hole 
doping level, the Bi-O band falls below Ep and, as a re- 
sult of this self-doping effect, further reduction in the hole 
doping level no longer reduces the number of holes in the 
Cu02 layers. We argue that the underlying mechanism 
at play here is that hole doping reduces the effective pos- 
itive charge in the Bi-O donor layers, which then reduces 
the tendency of the electrons to 'flow back' and self-dope 



the material. We have also carried out computations 
on a number of related compounds, including monolayer 
and trilayer Bi-compounds [Bi2Sr2Cu06+5 (Bi2201) and 
Bi2Sr 2 Ca2Cu30io +5 (Bi2223)], and the Tl- and Hg-based 
compounds, and we find that the lifting of the cation- 
derived band with hole doping is a generic property of 
many families of cuprates. 

Concerning technical details, we have employed both 
the Korringa-Kohn-Rostoker (KKR) and linearized aug- 
mented plane wave (LAPW) band structure method- 
ologies where we treat all electrons in the system self- 
consistently and consider the full crystal potential with- 
out the muffin-tin approximation^2*ii The KKR scheme 
is well-known to be particularly suited for a first- 
principles treatment of the electronic structure of sub- 
stitutionally disordered alloys. Pb substitution on the 
Bi sites was considered within the framework of the vir- 
tual crystal approximation (VCA), where the Bi nuclear 
charge Z is replaced by the average of the Bi and Pb 
charges of what may be thought of as an 'effective' Bi/Pb 
atom, but otherwise the band structure problem is solved 
fully self-consistently maintaining the charge neutrality 
of the system. The VCA is expected to be a good ap- 
proximation in this case since the effective disorder pa- 
rameter for the Bi-O states, given by A/W, where A is 
the splitting of the Bi-O and Pb-0 bands in Bi22I2 and 
Pb22I2, respectively, and W is the band width, is esti- 
mated to be ~ 0.3, so that the system is far from being 
in the split-band limit ^ We have also carried out su- 
perlattice computations by substituting two Bi atoms by 
Pb in the orthorhombic Bi22I2 as well as KKR-CPA (co- 
herent potential approximation) computations^ in 10% 
Pb doped Bi2201 to independently verify that the VCA 
provides a good description and that Bi/Pb substitution 
causes little disorder induced smearing of states. 

We set the stage for our discussion by considering 
Fig. 1(a) which shows the band structure of undoped 
Bi2212 predicted by the LDA. Here the lattice is as- 
sumed to be tetragonal and the structural parameters 
used are those of Ref. |U obtained by minimizing the 
total energy. A pair of closely placed bands is seen to 
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FIG. 1: Band structure of undoped and 25% Pb-doped 
tetragonal Bi2212 along various high symmetry directions at 
k z = 0. 



disperse rapidly through Ep along the T — X(ir, ir) line 
on the right side of Fig. 1(a). These are the well-known 
CuC>2-bands which are split into bonding and antibond- 
ing combinations due to intracell interactions between 
the two Cu02-planes. The problem however is that ad- 
ditional bands of BiO-character drop below Ep at the 
M(ir, 0) point giving the so-called 'Bi-0 pockets', lead- 
ing to a metallic Bi-0 layer ^ in clear disagreement with 
experimental observations^. 

Fig. 1(b) shows how the band structure changes dra- 
matically around the M-point when 25% Pb is substi- 
tuted for Bi in Bi2212, where the band structure of the 
doped compound is computed within the VCA^ The 
Bi-0 pocket problem is cured as the Bi-0 bands are 
lifted to Ri 0.4 eV above Ep, and the band structure 
around the M-point is simplified and the bilayer split- 
ting of the Cu02 bands around M becomes more clearly 
visible. The extended van Hove singularities (VHSs) in 
the antibonding and bonding bands appear at binding 
energies of -0.07 eV and -0.45 eV, respectively, and the 
bare bilayer splitting at M is « 400 meV. The shape of 
the antibonding and bonding Cu02-bands is very simi- 
lar to the generally accepted shape in the cuprates. The 
bands in Fig. 1(b) closely resemble the bands obtained in 
previous computations^*^ where an ad hoc modification 
of the LDA potential was invoked to account for the ab- 
sence of Bi-0 pockets in the ARPES spectra of Bi2212. 
Further computations for a range of Pb-doping levels in- 
dicate that the Bi-0 pockets are lifted just above Ep at 
around 22% Pb-doping in the tetragonal structure. 

The crystal structure of Bi2212 is more realistically 
modeled as a \pl x \/2 orthorhombic unit celliiV Accord- 
ingly, Fig. 2 delineates the effect of doping in orthorhom- 
bic Bi2212iSiiLi^ Here one obtains twice the number of 
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FIG. 2: Band structure (at fc z = 0) of undoped and 15% 
Pb-doped Bi2212 assuming orthorhombic lattice structure. 
Bands are plotted along the high symmetry lines of the tetrag- 
onal lattice for ease of comparison with the results of Fig. 1. 



bands compared to tetragonal Bi2212 due to the larger 
size of the unit cell. Comparing Figs. 2(a) and (b), we see 
that 15% Pb doping in the orthorhombic case lifts the Bi- 
O pockets « 0.4 eV above Ep and yields a FS consisting 
of only the bonding and antibonding Cu02 sheets. 

A comparison of Figs. 1 and 2 reveals interesting dif- 
ferences between the band structures of tetragonal and 
orthorhombic Bi2212 and their evolution with Pb doping. 
The Bi-0 complex of bands is more spread out in energy 
in Fig. 2(a) than in 1(a), which reflects the larger atomic 
displacements in Bi-0 layers in the orthorhombic struc- 
ture. The Bi-0 bands display greater sensitivity to Pb 
doping in the orthorhombic case and only 12% Pb doping 
pushes the Bi-0 pockets above Ep compared to the value 
of 22% needed in the tetragonal structure. There also are 
differences in the Cu02 bands. For example, the doping 
level at which the VHS of the antibonding band lies at 
the Ep is 22% in the orthorhombic structure and 27% in 
the tetragonal case. Besides the highly dispersive Cu02 
bands, the complex of Cu-0 bands below Ep (starting 
around a binding energy of sa 0.8 eV in Fig. 2(a)), which 
is primarily composed of Cu d and O p bands, is also 
influenced by the crystal structure and doping as seen 
with reference to Figs. 1 and 2. In particular, in the 
doped orthorhombic system in Fig. 2(b), around the M- 
point, these lower lying bands mix significantly with the 
Cu02 band involved in producing the bonding FS sheet, 
and change the shape of the associated VHS. 

Fig. 3 shows that our theoretically predicted FS is in 
remarkable accord with the experimentally determined 
FS of an overdoped Bi2212 sample obtained via angle- 
resolved photoemission (ARPES) measurements 9 . For 
this purpose, computed FS contours for Pb doping levels 
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FIG. 3: (Color) Theoretical bonding (yellow lines) and an- 
tibonding (red lines) FS's of orthorhombic Bi2212 for 20% 
(solid lines) and 30% (dashed lines) Pb doping. Other 
'shadow' FS's are not shown for simplicity. Experimental FS 
map taken via ARPES from an overdoped single crystal of 
Bi2212 is from Ref.@. 



of 20% (solid lines) and 30% (dashed lines) for the or- 
thorhombic lattice are overlayed on the experimental FS 
map. 19,20,21 The 'shadow' FS's in the computations arc 
not shown in order to highlight the main bonding and an- 
tibonding FS's. The computed bonding FS (yellow lines) 
shows relatively little change over 20-30% doping range 
and its shape and dimensions are in quantitative accord 
with measurements. The antibonding FS (red lines), on 
the other hand, is more sensitive to doping and changes 
from being hole-like at 20% doping (solid red line) to 
turning electron-like (dashed red line) at 30% doping as 
the Ep descends through the VHS. Therefore the spec- 
tral intensity associated with the antibonding FS in the 
antinodal region will be sensitive to local variations in 
hole doping and a careful modeling of the spectral in- 
tensities will be required to pin down details of the FS. 
However, along the nodal direction, neither the antibond- 
ing nor the bonding FS is sensitive to doping and here 
there is good accord between theory and experiment. 

The driving mechanism underlying the lifting of the Bi- 
O pockets with Pb doping in our computations may be 
understood as follows. The ionization of Bi atoms in the 
system will in general generate electric fields which tend 
to attract electrons back into the Bi-0 layers and com- 
pete with the affinity for the electrons towards the CuC>2 
layers. The band structures of Figs. 1(a) or 2(a) which 
display partially filled Bi-0 bands and the associated Bi- 
O pockets at the FS then imply that the balance of forces 
in the computations is such that Bi is not fully ionized to 
3+ in pure Bi2212 so that we may think of some of the 
Bi 3+ electrons as being attracted back to the Bi-0 layers 



FIG. 4: Band structures of O-doped Bi2212: (a) S = 0.1, and 
(b) S = 0.3, where S denotes excess O per unit cell. 



or that the CuC>2 layers are self-doped with holes. The 
fact that the Bi-0 bands are moved above Ep with Pb 
doping in Figs. 1(b) or 2(b) then indicates that the sub- 
stitution of Bi with Pb and the concomitant reduction 
of positive charge in the Bi-0 layers eliminates the need 
for electrons to 'flow back' to the Bi-0 layer. In effect 
then at e.g. 25% Pb doping an empty (Bi/Pb)-0 band 
only donates 0.75 rather than 1.0 electron to the Cu02 
layer. It is helpful as well to see how this argument plays 
out in reverse, i.e. with decreasing Pb doping. When Pb 
doping decreases, the tendency of the Bi/Pb electrons to 
flow back to the (Bi/Pb)-0 layer increases and below a 
critical Pb-doping level some of the Bi/Pb electrons ac- 
tually flow back as the (Bi/Pb)-0 band drops below Ep. 
As a result of this self-doping effect, further decrease in 
hole doping of the Cu02 layers is prevented. 

We turn now to consider the effect of adding excess 
O in Pb-free Bi2212 for hole doping the system. For 
this purpose, we have carried out extensive computations 
where O, F, or other pseudo-atoms are inserted in the 
empty spaces between the Bi-0 layers in order to cap- 
ture varying amounts of Bi electrons to form a closed 
shell 2 ^. Typical modifications in the band structure are 
shown in Fig. 4, where we see that the effect of excess 
O is to lift the Bi-0 pockets much like that of Pb/Bi 
substitution. The key is to reduce the effective number 
of electrons available in the Bi-0 layers for donation to 
the Cu02 layers and this can be accomplished via cither 
Pb/Bi substitution or by adding excess O^^Sk In Fig. 
4, the Bi-0 pockets lie below Ep for excess O value of 
<5= 0.1, but lie well above Ep for 6— 0.3. Our analysis 
indicates that the Bi-0 pockets move through the Ep at 
5w 0.18. 

We have also considered the effect of hole doping on 
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other Bi-compounds and find that the Bi-0 bands in the 
15% Pb doped Bi2201 and Bi2223 are lifted above E F , 
even though the band structures of the undoped com- 
pounds in both cases display Bi-0 FS pocketsi^Si2L2& 
Going beyond the Bi-compounds, we have studied dop- 
ing effects on the band structures and binding energies 
of the core levelsS^H in the Tl- and Hg-based cuprates. 
In particular, the Tl-0 pockets around T-point in Tl- 
based cuprates are removed by approximately 10% hole 
doping^i Interestingly, in Tl 2 Ba 2 Cu0 6 (T12201), in- 
clusion of hole doping results in a 3D FS that agrees 
well with angular magnetoresistance oscillation measure- 
ments of Ref. |32j The Hg-derived bands in the Hg-based 
cuprates are similarly lifted by O-doping. 

In order to assess the effect of superlattice modula- 
tion, we have carried out a series of computations in the 
orthorhombic structure where the refined parameters of 
Miles et ali& at various positions along the modulation 
were used. We found that the Bi-0 pockets could not 
be lifted for any value of the parameters'^ Even oth- 
erwise, we would not expect superlattice modulation to 
provide a generic mechanism for lifting the cation-derived 
bands in the cuprates because these modulations vary 
greatly between different cuprates, and are suppressed 
by Pb doping in the Bi-compounds. 



In conclusion, our results show clearly that substantial 
and generic Coulombic effects come into play with hole 
doping to lift the cation-derived bands in the cuprates. In 
adducing various physical quantities from spectroscopic 
data (e.g. size of the pseudogap), changes in the elec- 
tronic structure with underdoping, especially near the 
anti-nodal point, should be accounted for, even though 
most of the existing analysis in the cuprate literature as- 
sumes a doping-independent band structure. Finally, our 
study provides a first-principles route for exploring self- 
doping effects and doping dependencies of the electronic 
structures of this exciting class of materials. 
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